perm filename DEG6.IN[225,JMC] blob
sn#005373 filedate 1971-03-05 generic text, type T, neo UTF8
00100 MATRIX MZ(5,5),M12(5,5),M34(5,5),M23(5,5),M45(5,5),M56(5,5);
00200
00300 FOR I←1 STEP 1 UNTIL 5 DO
00400 FOR J←1 STEP 1 UNTIL 5 DO
00500 MZ(I,J)←0;
00600
00700 M12←M34←M23←M45←M56←MZ;
00800
00900 M56(1,1)←C1;M56(2,2)←C2;M56(3,3)←C3;M56(4,4)←C4;M56(5,5)←C5;
01000
01100 M45←MAT((X11,X12,X13,X14,X15),(X21,X22,X23,X24,X25),
01200 (X31,X32,X33,X34,X35),(X41,X42,X43,X44,X45),(X51,X52,X53,X54,X55));
01300
01400 M12(4,4)←M12(5,5)←M23(4,4)←M23(5,5)←M34(4,4)←M34(5,5)←1;
01500
01600 FOR I←1:3 DO FOR J←1:3 DO BEGIN
01700 M12(I,J)←B12(I,J);M23(I,J)←B23(I,J);M34(I,J)←B34(I,J) END;
01800
01900 OFF ECHO;
02000
02100 END;
02200
02300 C1←C2←C3←C4←1; C5←-1;
02400
02500 M45←MAT((3/4,-1/2,0,-3/4,1),(-1/4,1/2,0,-3/4,1),
02600 (0,0,1,0,0),(-1/4,-1/2,0,1/4,1),(1/8,1/4,0,3/8,1/2));